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ABSTRACT 

We present ground-based optical observations of GRB 020124 starting 1.6 hours after the burst, as well 
as subsequent Very Large Array (VLA) and Hubble Space Telescope (HST) observations. The optical after- 
glow of GRB 020124 is one of the faintest afterglows detected to date, and it exhibits a relatively rapid decay, 
r-i *'0±o.04^ followed by further steepening. In addition, a weak radio source was found coincident with the 
optical afterglow. The HST observations reveal that a positionally coincident host galaxy must be the faintest 
host to date, R ^ 29.5 mag. The afterglow observations can be explained by several models requiring little or 

no extinction within the host galaxy, Ajj,™' w 0-0.9 mag. These observations have significant implications for 
the statistics of the so-called dark bursts (bursts for which no optical afterglow is detected), which are usually 
attributed to dust extinction within the host galaxy. The faintness and relatively rapid decay of the afterglow of 
GRB 020124, combined with the low inferred extinction indicate that some dark bursts are intrinsically dim and 
not dust obscured. Thus, the diversity in the underlying properties of optical afterglows must be observationally 
determined before substantive inferences can be drawn from the statistics of dark bursts. 
Subject headings: gamma-rays:bursts — dust:extinction — cosmology:observations 



1. INTRODUCTION 

One of the main observational results stemming from five 
years of 7-ray burst (GRB) follow-ups at optical wavelengths 
is that about 60% of well-localized GRBs lack a detected op- 
tical afterglow, ("dark bursts"; Taylor et al. 2000; Fynbo et 
al. 2001; Reichart & Yost 2001; Lazzati, Covino, & Ghis- 
eUini 2002). In some cases, a non-detection of the optical after- 
glow could simply be due to a failure to image quickly and/or 
deeply enough. However, there are two GRBs for which there 
is strong evidence that the optical emission should have been 
detected, based on an extrapolation of the radio and X-ray emis- 
sion (Djorgovski et al. 2001a; Piro et al. 2002). One interpreta- 
tion in these two cases is that the optical light was extinguished 
by dust, either within the immediate environment of the burst 
or elsewhere along the line of sight (e.g. Groot et al. 1998). An 



alternative explanation is a high redshift, leading to absorption 
of the optical light in the Lya forest. However, the redshifts of 
the underlying host galaxies of these GRBs are of order unity 
(Djorgovski et al. 2001a; Piro et al. 2002). 

Several authors have recently argued that a large fraction of 
the dark bursts are due to dust extinction within the local en- 
vironment of the bursts (e.g. Reichart & Yost 2001; Lazzati et 
al. 2002; Reichart & Price 2002), but other scenarios have also 
been suggested (e.g. Lazzati et al. 2002). Moreover, it has 
been noted that regardless of the location of extinction within 
the host galaxy, the fraction of dark bursts is a useful upper 
limit on the fraction of obscured star formation (Kulkarni et al. 
2000; Djorgovski et al. 2001b; Ramirez-Ruiz, Trentham, & 
Blain 2002; Reichart & Price 2002). 

However, from an observational point of view, we must have 
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a clear understanding of the diversity of afterglow properties be- 
fore extracting astrophy sic ally interesting inferences from dark 
bursts. For example, afterglows which are faint or fade rapidly 
(relative to the detected population) would certainly bias the 
determination of the fraction of truly obscured bursts. In this 
vein, Fynbo et al. (2001), noting the faint optical afterglow of 
GRB 000630, argued that some dark bursts are due to a failure 
to image deeply and/or quickly enough, rather than dust extinc- 
tion. 

Here we present optical and radio observations of 
GRB 020124, an afterglow that would have been classified dark 
had it not been for rapid and deep searches. Furthermore, 
GRB 020124 is an example of an afterglow, which is dim due 
to the combination of intrinsic faintness and a relatively fast 
decline, and not strong extinction. 

2. OBSERVATIONS 
2.1. Ground-Based Observations 

GRB 020124, locaHzed by the HETE-II satellite on 2002, Jan 
24.44531 UT, had a duration of ~ 70 s and a fluence (6-400 
keV) of 3 X 10"^ erg cm'^ (Ricker et al. 2002). Eight min- 
utes after receiving the coordinates'^ we observed the error 
box with the dual-band {Bm, Rm) MACHO imager mounted on 
the robotic 50-in telescope at the Mount Stromlo Observatory 
(MSO). We also observed the error box with the Wide-Field Im- 
ager on the 40-in telescope at Siding Spring Observatory (SSO). 
We were unable to identify a transient source within the large 
error box (Price, Schmidt, & Axelrod 2002). 

We subsequently observed the error box refined by the Inter- 
Planetary Network (Hurley et al. 2002) with the Palomar 48- 
in Oschin Schmidt using the unfiltered NEAT imager. PSF- 
matched image subtraction (Alard 2000) between the MA- 
CHO and NEAT images revealed a fading source (Price et al. 
2002), which was w 18 mag at the epoch of our first ob- 
servations, and not present in the Digitized Sky Survey. Two 
nights later we observed the afterglow using the Jacobs CAM- 
era (JCAM; Bloom et al. 2002b) mounted at the East arm 
focus of the Palomar 200-in telescope (Bloom et al. 2002a). 
The position of the fading source is a(J2000)=9''32™50.78^ 
(5(J2000)=-11°31'10.6", with an uncertainty of about 0.4 arc- 
sec in each coordinate (Fig. |l]). 

Using the Very Large Array (VLA'^) we observed the fad- 
ing source at 8.46 and 22.5 GHz (see Table |]). We de- 
tect a faint source, possibly fading, at 8.46 GHz located at 
a(J2000)=9^32"50.8P, (5(J2000)= -11°31'10.6", with an un- 
certainty of about 0.1 arcsec in each coordinate. Given the 
positional coincidence between the fading optical source and 
radio detection we suggest this source to be the afterglow of 
GRB 020124. 

The optical images were bias-subtracted and flat-fielded in 
the standard manner To extract the photometry we weighted 
the aperture with a Gaussian equivalent to the seeing disk 
("weighted-aperture photometry"), using IRAF/wphot. The 
photometric zero-points were set through photometry of cali- 
brated field stars (Henden 2002) with magnitudes transformed 
to the appropriate system (Bessell & Germany 1999; Smith et 
al. 2002). The photometry is summarized in Tab. [l| 

"This corresponds to 1.6 hours after the burst detection. 
'*The VLA is operated by the National Radio Astronomy Observatory, a fi 
Associated Universities, Inc. 



2.2. Hubble Space Telescope Observations 

We observed the afterglow with the Hubble Space Telescope 
(HST) using the Space Telescope Imaging Spectrograph (STIS) 
on 2002 Feb. 11.09, 18.30, and 25.71 UT (Bloom et al. 2002a), 
as part of our large HST Cycle 10 program (GO-9180, PI: 
KuUcarni). The HST observations consisted of 750-850 sec ex- 
posures. The HST data were retrieved after "On-The-Fly" pre- 
processing. Using IRAF we drizzled (Fruchter & Hook 2002) 
each image onto a grid with pixels smaller than the original by 
a factor of two and using pixf rac of 0.7. 

We found an astrometric tie between the HST and JCAM 
images using IRAF/geomap with nine suitable astrometric tie 
objects in common between the images. The rms of the resul- 
tant mapping is 133 mas (RA) and 124 mas (Dec). Using this 
mapping and IRAF/geoxytrans we transfered the afterglow 
position on the JCAM image to the HST images. The rms of 
the transformation is 604 mas (RA) and 596 mas (Dec), and is 
dominated by the uncertainty in the JCAM position. 

The source "SI" (Fig. ^ coincides with the afterglow po- 
sition within the astrometric uncertainty. We performed dif- 
ferential photometry at the position of SI by registering the 
images of epochs 1 and 2 using a cross-correlation of a field 
of size 10 arcsec centered on SI (using IRAF/crosscor and 
shift find). We used IRAF/center and the FWHM of a 
relatively bright point source ("PSF star"; Fig. |l]) to fix the po- 
sition of SI in each of the final images, and to determine the 
uncertainty in the position. 

We photometered the source (and the PSF star) in epoch 1 us- 
ing IRAF/phot, and determined a count-rate of 0.0921 ±0.013 
e" s"' (corrected by 17% for the loss of flux from an infinite 
aperture radius). Using IRAF/synphot and assuming a source 
spectrum of f\ oc A"' '* (see below), we find that the source was 
R = 28.55;!;o [4 '■i™^ °f epoch 1. The photometry of 

the three epochs is summarized in Tab. |[ Please note that this 
more careful analysis supersedes our preliminary report (Bloom 
et al. 2002c). 

There are no obvious persistent sources within 1.75 arcsec of 
the OT down to w 29.5 mag. To date, aU of the GRBs local- 
ized to sub-arcsecond accuracy have viable hosts brighter than 
this level within 1 .3 arcsec of the OT position (Bloom, KuUca- 
rni, & Djorgovski 2002). The faintest host to date is that of 
GRB 990510, R - 28.5 mag (z=1.619; Vreeswijk et al. 2001). 
Thus, the host of GRB 020124 may be at a somewhat higher 
redshift; however, z <~ 4.5 since the afterglow was detected in 
the Bm filter. 

3. MODELING OF THE OPTICAL DATA 

In Figure ^ we plot the optical lightcurves of GRB 020124, 
including a correction for Galactic extinction, E{B-V) = 0.052 
mag (Schlegel et al. 1998). The optical lightcurves are usu- 
ally modeled as F^{t,v) = Fy x){t / tQ)°'{v / vq)^ . However, as can 
be seen in Fig. ^ the /?-band lightcurve cannot be described 
by a single power law. Restricting the fit to f < 2 days we ob- 
tain (Xmin - 15 foJ" 14 degrees of freedom) a\ =-1.60 ±0.04, 
13 = -1.43 ±0.14, and /v.o = 2.96 ±0.25 /^Jy; here /v^) is de- 
fined at the effective frequency of the Rm filter and f = 1 day. 
For f > 2 days we get 02 « -1 .9. 

of the National Science Foundation operated under cooperative agreement by 
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To account for the steepening we modify the model for the 
/?-band lightcurve to: 



We find fjet ^ 10-20 days, corresponding to 6'jet ^ 10°. Using 



the measured fluence (§ 2. 1 ) we estimate the beaming-corrected 



F,(f , y) = F,.o(y/iyof[(t/tbr" + (t/thm 



l/n 



(1) 



where, ai is the asymptotic index for t <C ?/,, 02 is the asymp- 
totic index for f » f^,, n < provides a smooth joining of the 
two asymptotic segments, and tb is the time at which the asymp- 
totic segments intersect. We retain the simple model for the Rm 
and Bm lightcurves since they are restricted to f < 0.13 days 
(i.e. well before the observed steepening). 

We investigate two alternatives for the observed steepening 
in the framework of the afterglow synchrotron model (e.g. Sari, 
Piran, & Narayan 1998). In this framework, ai, a2, and f3 are 
related to each other through the index (p) of the electron en- 
ergy distribution, A^(7) c>c 7"'' (for 7 > 7min). The relations for 
the models discussed below, as well as the resulting closure re- 
lations, ai+b/S + c = 0, are summarized in Tab. ^ 

3.1. Cooling Break 

The observed steepening, Aa = 0,2- oi\ « -0.3, can be due 
to the passage of the synchrotron cooling frequency, Vc, through 
the /?-band. This has been inferred, for example, in the after- 
glow of GRB 971214, at t - 0.2 days (Wijers & Galama 1999). 
If the steepening is due to Vc, this rules out models in which 
the ejecta expand into a circumburst medium with p oc r"^ 
(hereafter. Wind), because in this model Vc increases with time 
(oc f Chevalier & Li 1999), and one expects Aa = 0.25. 

There are two remaining models to consider in this case: (i) 
spherical expansion into a circumburst medium with constant 
density (hereafter, ISMb; Sari, Piran, & Narayan 1998), and 
(ii) a jet with 6'jet < r7~0 06d (i-^- ^ J^*- br^^k prior to the first 
observation at f « 0.06 days; hereafter. Jets). The subscript B 
indicates that is blueward of the optical bands initially. In 
both models we use Eqn. |l] for the /?-band lightcurve, with f/, 
defined as the time at which Vc = vr, and a2 = a\-\ /4. 

We find that in the ISMb model tc « 0.4 days, while in the 
Jetfi model tc ~ 0.65 days. Moreover, in both models the clo- 
sure relations can only be satisfied by including a contribution 
from dust extinction within the host galaxy. Ay"'*'. We esti- 
mate the required extinction using the parametric extinction 
curves of Cardelli, Clayton, & Mathis (1989) and Fitzpatrick 
& Massa (1988), along with the interpolation calculated by Re- 
ichart (2001). Since the redshift of GRB 020124 is not known 
we assume z = 0.3, 1, 3, which spans the range of typical red- 
shifts for the long-duration GRBs. The inferred values of A^^^ 
are summarized in Tab. ^ and range from 0.2 to 0.9 mag. 

3.2. Jet Break 

An alternative explanation for the steepening is a jet expand- 
ing into: (i) an ISM medium with Vc blueward of the optical 
bands (J-ISMb), (ii) a Wind medium with Vc blueward of the 
optical bands (J- Winds), and (iii) an ISM or Wind medium with 
Vc redward of the optical bands (J-ISM/WindR). We note that 
the J-ISMb model is different than the ISMb model (§01) since 



previously it was defined such that the jet break is later than the 
last observation. In these models, tt, = fjet is the time at which 

From the closure values we note that the J-ISM/WindR re- 
quires no extinction within the host galaxy, while the J-ISMb 
and J-Windfi models require values of about 0.05 to 0.3 mag. 



7-ray energy, E-^ ^5 x 10 erg, assuming a circumburst den- 
sity of 1 cm"-' and z = 1 {E,y is a weak function of z). This 
value is in good agreement with the distribution of E,y for long- 
duration GRBs (Frail et al. 2001). 

4. DISCUSSION AND CONCLUSIONS 

Regardless of the specific model for the afterglow emis- 
sion, the main conclusion of §|^ is that the optical afterglow of 
GRB 020124 suffered litde or no dust extinction. Still, this af- 
terglow would have been missed by typical searches undertaken 
even as early as 12 hours after the GRB event. As shown in 
Fig.0 about 70% of the searches conducted to date would have 
failed to detect an optical afterglow like that of GRB 020124. 

This is simply because the afterglow of GRB 020124 was 
faint and exhibited relatively rapid decay. From Fig. ^ we note 
that GRB 020124 is one of the faintest afterglows detected to 
date (normalized to f = 1 day), and while it is not an excessively 
rapid fader, it is in the top 30% in this category. 

Thus, the afterglow of GRB 020124^ along with that of 
GRB 000630 (Fynbo et al. 2001; Fig. ||), indicates that there 
is a wide diversity in the brightness and decay rates of optical 
afterglows. In fact, the brightness distribution spans a factor of 
about 400, while the decay index varies by more than a factor of 
three. Coupled with the low dust extinction in the afterglow of 
GRB 020124, this indicates that some dark bursts may simply 
be dim, and not dust obscured. 

Given this wide diversity in the brightness of optical after- 
glows, it is important to establish directly that an afterglow 
is dust obscured. This has only been done in a few cases 
(§|ll)- Therefore, while statistical analyses (e.g. Reichart & Yost 
2001) point to extinction as the underlying reason for some 
fraction of dark bursts, and may even account for an afterglow 
like that of GRB 020124, it is clear that observationally the is- 
sue of dark bursts is not settled, and the observational biases 
have not been traced fully. 

Since progress in our understanding of dark bursts will ben- 
efit from observations, we need consistent, rapid follow-up of 
a large number of bursts to constrain the underlying distribu- 
tion, as well as complementary techniques which can directly 
measure material along the line of sight. This includes X-ray 
observations which allow us to measure the column density to 
the burst (Galama & Wijers 2001), and thus infer the type of 
environment, and potential extinction level. Along the same 
line, radio observations allow us to infer the synchrotron self- 
absorption frequency, which is sensitive to the ambient density 
(e.g. Sari & Esin 2001); the detection of radio emission, as in 
the case of GRB 020124, implies a density n <^ 10^ cm"^. Fi- 
nally, prompt optical observations, as we have carried out in 
this case, may uncover a larger fraction of the dim optical af- 
terglows, and provide a better constraint on the fraction of truly 
obscured bursts. 
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Table 1 

Ground-Based Optical Observations of GRB 020124 



UT 


Telescope 


Band 


Magnitude 


Jan 24.51204 


MSO50 


Rm 


17.918 ±0.041 


Jan 24.51204 


MSO50 


Bm 


18.628 ±0.057 


Jan 24.51516 


SSO40 


R 


18.219 ±0.046 


Jan 24.51655 


MSO50 


Rm 


17.984 ±0.044 


Jan 24.51655 


MSO50 


Bm 


18.727 ±0.063 


Jan 24.51938 


SSO40 


R 


18.371 ±0.091 


Jan 24.52106 


MSG 50 


Rm 


18.111 ±0.049 


Jan 24.52106 


MSG 50 


Bm 


18.842 ±0.069 


Jan 24.52373 


SSO40 


R 


18.376 ±0.082 


Jan 24.55791 


MSG 50 


Rm 


18.678 ±0.048 


Jan 24.55791 


MSG 50 


Bm 


19.661 ±0.090 


Jan 24.56243 


MSG 50 


Rm 


18.867 ±0.036 


Jan 24.56243 


MSG 50 


Bm 


19.584 ±0.053 


Jan 24.56696 


MSG 50 


Rm 


18.843 ±0.039 


Jan 24.56696 


MSG 50 


Bm 


19.7 14 ±0.050 


Jan 26.34100 


P200 


r' 


24.398 ± 0.228 



Note. — The columns are (left to right), (1) UT date of each observation, (2) telescope (MSO 50: Mt. Stromlo Observatory 50-in; SSO 40; Siding Spring Obser- 
vatory 40 in; P 200: Palomar Observatory 200-in), (3) observing band, and (4) magnitudes and uncertainties. The observed magnitudes are not corrected for Galactic 
extinction. 
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Table 2 

HST/STIS Observations of GRB 020124 



Epoch 
(UT) 



Band 



Exp. Time 
(ksec) 



Flux 

(e-s-l) 



S/N 



Magnitude 



Feb 11.09 
Feb 18.30 
Feb 25.71 



50CCD/Clear 
50 CCD/Clear 
50CCD/Clear 



10.0 
7.4 
7.5 



0.0814±0.0169 
0.0443 ±0.0189 
0.0362 ±0.0183 



4.82 
2.34 
1.98 



R — 28.68^„~„ 
i? = 29.35!g'f 
« = 29.56||. 



Feb 18.30+25.71 



50CCD/Clear 



14.9 



0.0398 ±0.0137 



2.91 



R = 29.46 



Note. — The columns are (left to right). (1) UT date of each observation. (2) STIS CCD mode, (3) exposure time, (4) flux and uncertainty significance, and (6) 
R magnitude and uncertainty. The total number of counts was converted to the R-band assuming the observed color of the OT, fx cx; A"" (§E.2|). For epochs 2 and 3, 
the 3(T upper limits are: R = 29.09 and R = 29.13 mag, respectively. The observed magnitudes are not corrected for Galactic extinction. 
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Table 3 

VLA Radio Observations of GRB 020124 



Epoch 




Flux Density 


(UT) 


(GHz) 


(MJy) 


Jan 26.22 


8.46 


84 ±30 


Jan 26.25 


22.5 


-60 ±100 


Jan 27.22 


8.46 


45 ±25 


Feb 1.40 


8.46 


49 ±17 


Jan 26.22-Feb 1.40 


8.46 


48 ±13 



Note. — The columns are (left to right), (1) UT date of each observation, (2) observing frequency, and (3) flux density at the position of the radio transient with the 
rms noise calculated from each image. The last row gives the flux density at 8.46 GHz from the co-added map. 
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Table 4 
Afterglow Models 



Model 


ai 


012 


P 


(b,c) 


Closure 


P 


A*™' (mag) 


ISMb 


3(p-l) 
4 


4^2 




(-3/2,0) 


0.52 ±0.28 


3.17±0.05 


(0.35,0.18,0.10) 


JetB 


-P 


-P 




(-2,1) 


2.23 ±0.36 


1.63 ±0.04 


(0.89, 0.50, 0.22) 


J-1SMb 


■i(p-l) 


-P 


A 


(-3/2,0) 


0.52 ±0.28 


3. 17 ±0.05 


(0.30, 0.10, 0.05) 


J-WindB 




-P 


2 


(-3/2,1/2) 


1.02 ±0.28 


2.51 ±0.05 


(0.30, 0.16, 0.08) 


J-lSM/WindR 


3p-2 
4 


-P 


P 
2 


(-3/2,-1/2) 


0.02 ±0.28 


2.84 ±0.05 





Note. — The columns are (left to right), (1) Afterglow model (ISM: r" circumburst medium; Wind: r^^ circumburst medium; Jet: collimated eject with opening 
angle 6jet; a subscript B indicates i/c < i^opt, and a subscript R indicates > fopt), (2) aj as a function of p, (3) a2 as a function of p, (4) /3 as a function of/), (5) 
closure relations (a + bp + c = 0), (6) resulting closure values from the observed values of a[ and /3, (7) inferred value of p from the measured value of cti, and (8) 
the required extinction in the frame of the host galaxy for closure values of zero (z = 0.3, 1, 3); typical uncertainties are ±0.05 mag. The top two models apply to the 
case when the observed steepening in the lightcurves is due to the passage of Uc through the R-band, while the bottom three apply to the case when the steepening is 
due to a jet. 
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Fig. 1. — Palomar 200-inch (left) and HST epoch 1 (inset) images of the field of GRB 020124. The OT is circled in both images. 
The OT was of comparable brightness to Gl at the epoch of the P 200 image and significantly fainter than Gl three weeks later The 
box overlaying the inset shows the portion of the HST images depicted in Figure ^ Relevant sources described in the text are noted. 
The HST image is shown with logarithmic scaling to highlight the features of nearby galaxies.. 
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Fig. 2.— The faint optical transient (OT) of GRB 020124 as viewed using HST/STIS. Shown are the summed, smoothed images 
from epoch 1 (left) and epochs 2+3 (right). The greyscales have been matched such that a given flux is represented by the same shade 
in each image. The circle is centered at the same sky position in both imagess. Clearly, the source SI, identified with the position of 
the afterglow of GRB 020124 has faded. 
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Time since the burst (days) 

Fig. 3.— Optical lightcurves of GRB 020124 (top to bottom: Rm, R, and Bm), corrected for Galactic extinction, E{B-V) = 0.052 
mag (Schlegel et al. 1998). The solid lines are a representative jet model (ISM/Wind^; see Tab. while the dashed line is an 
extrapolation of the early evolution without a break. With no break in the /?-band lightcurve, the predicted magnitude at the epoch 
of the first HST observation exceeds the measured values by 5a. The flux measured in the last HST epoch is plotted as a 2t7 upper 
limit. 



12 



18 


GRB 020124 


GRB 970828 \gRB 991208 
\ (pred. from X-rays) ' \ 


- 


20 


GRB 000630 


\. V V 




cc 
cc 




\. V 

V ^^^^^^ \. V \ ^ 

\^ V \ 




22 










GRB 970828 ' ^^"^f^i 
(obs. limit) w 




24 


1 


V V 
V \^ 





0.1 1 

Time After the Burst (Days) 



Fig. 4. — R-hand upper limits from searches of well-localized GRBs, corrected for Galactic extinction. The limits up to GRB 000630 
are taken from Fynbo et al. (2001), while subsequent limits are from the GRB Coordinates Network. Also shown are the lightcurves 
of the GRB 020124, GRB 000630, the bright GRB 991208 (Castro-Tirado et al. 2001), and GRB 970828 (the de-reddened Ughtcurve 
is based on the radio and X-ray data; Djorgovski et al. 2001). Only about 30% of the searches yielded limits that are fainter than the 
afterglow of GRB 020124. A similar fraction was found by Fynbo et al. (2001) based on the afterglow of GRB 000630. 
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Fig. 5. — Temporal decay index, ckopt {Fu oc f"), plotted against the R-hwA magnitude at f = 1 day for several optical afterglows. 
We chose a fiducial time of 1 day since, with the exception of GRB 010222, all the observations are before the jet break. While the 
majority of optical afterglows cluster around R{t = 1 d) ~ 20 mag, GRB 020124 is one of the four faintest afterglows detected to date, 
and one of the six most rapid faders. 



